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Abstract. Guided by the Hopf fibration, we single out a family 
(indexed by a positive constant K) of right invariant Riemannian 
metrics on the Lie group S 3 . Using the Yasuda-Shimada theorem 
as an inspiration, we determine for each K > 1 a privileged right 
invariant Killing field of constant length. Each such Riemannian 
metric pairs with the corresponding Killing field to produce a y- 
global and explicit Randers metric on S 3 . Using the machinery of 
spray curvature and Berwald's formula for it, we prove directly that 
the said Randers metric has constant positive flag curvature K, as 
predicted by the Yasuda-Shimada theorem. We also explain why 
this family of Finslerian space forms is not projectively fiat. 



1. Introduction 

A Finsler metric F is a family of 'norms' on a manifold M, one on each 
tangent space T X M. These 'norms' are typically only positively homoge- 
neous of degree one, whereas the norms used in functional analysis are 
absolutely homogeneous. There are also the usual smoothness and strong 



convexity assumptions (see for instance BCS]) on the slit tangent bundle 
TM \ 0. In a large number of examples, especially the ones of physi- 
cal origin, these technical requirements are only satisfied on open cones in 
TM \ 0. If the Finsler metric is smooth and strongly convex on the entire 
slit tangent bundle TM \ 0, it is said to be y-global. Examples that are 
both y-global and geometrically significant are highly sought after. 

In Riemannian geometry, one has the concept of sectional curvature. 
Its analogue in Finsler geometry is called the flag curvature. Flag curva- 
tures are more easily accessed in some settings through spray curvatures. 
These objects have now been gi ven d etailed treatments in textbooks and 
monographs such as §§, §3CS|, lAIMfl , and Finsler spaces of constant 



flag curvature are, just like their Riemannian counterparts, known as space 
forms. However, unlike Riemannian geometry, Finslerian space forms do 
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not just arise from three standard models. This phenomenon is related to 
the fact that on K™, there are many isometry classes of norms, but only 
one isometry class of inner products. 

Akbar-Zadeh |AZ| showed that if the metrics in question are geodesically 
complete and the growth of their Cartan tensors are suitably constrained, 
then flat and negatively curved space forms are fairly well understood. See 
[BCS| for a leisurely exposition. Relax any of those two hypotheses and 
one encounters intriguing spaces, like for example the Finslerian Poincarc 



disc discussed in [BCS 



This paper deals with y-global Finsler spaces of constant positive flag 
curvature, and puts us in an even more esoteric landscape. In this regard, 
we have several pioneering works of Bryant's. Two dimensional examples 
are treated in [Brl], [Br2|, [Br3|; higher dimensional ones are discussed in 
his lectures. As far as we know, all of Bryant's examples are projectively 
flat (meaning that their geodesies are straight lines in certain coordinate 
systems), and none of them is of the Randers vintage. 

Randers spaces are Finsler spaces constructed from just two pieces of 
familiar data: a Ricmannian metric and a differential 1-form, both globally 
defined on an underlying smooth manifold. As such, they are possibly the 
best stepping stones from the Riemannian realm to the Finslerian territory. 

The goal of this paper is to produce, for each constant K > 1, an 
explicit example of a compact boundaryless (non-Riemannian) 
Randers space that has constant positive flag curvature K, and 
which is not projectively flat. 

For that purpose, we turn to the Yasuda-Shimada theorem Y£]. This 
result was published in the late 70s, and classifies Randers metrics of con- 
stant positive and constant negative flag curvature. Incidentally, flat Ran- 
ders metrics are necessarily locally Minkowskian. For the positively curved 
case, the Yasuda-Shimada theorem gives four mathematical criteria on the 
Riemannian metric and the 1-form. These criteria are made precise in §2. 

For each K > 1, we select a manifold M judiciously, then solve 
those four criteria for a Riemannian metric a and a nonzero 
1-form b, both living on M. This is carried out in §3— §5. By 
doing so, we have shown that the Yasuda-Shimada criteria are 
non-vacuous. The simplest choice of M turns out to be the Lie 
group S 3 . Similar promise [Ro] holds for 5 2ll+1 , mainly because 
every odd-dimensional sphere admits the Hopf fibration. 

After the construction, we could in principle appeal to the Yasuda- 
Shimada theorem to conclude that our Randers spaces have constant pos- 
itive flag curvature. However, the two published proofs ( JYSH and (mJ) of 
the theorem in question both contain arguments that we find difficult to 
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follow. In order to render our paper logically self-contained, we prove di- 
rectly that each said example indeed has constant positive flag curvature 
K. The method we use is quite different from the ones in 

flU and ||. 

It is geometrical and brings out the spectacular power behind a formula of 
Berwald's. It also reaffirms our belief that a happy synergy awaits Finsler 
geometry and modern computing. 

To that end, we first review the concept of spray curvatures, and its 
close relationship with the flag curvature, in §6. We then give Berwald's 
formula for the spray curvature. This formula has proved to be exception- 
ally useful in computational Finsler geometry, whether one carries out the 
computations on machine or by hand. 

In the Appendix, we write a Maple program based on Berwald's formula 
and use it to generate numerical data, supporting the contention that our 
Randers metrics have constant positive flag curvature K. This type of 
experimental evidence provides the faith which sustains our direct proof 
(by hand) in §7. The hand computations in this direct proof are facilitated 
by a refined understanding [ BjJ of the spray curvature, when the latter is 
adapted to the Randers setting. 

We believe that this refined understanding also holds the key to a dis- 
tinctly new proof of the Yasuda-Shimada theorem. This theorem has been a 
reliable source of beautiful examples. For instance, the Finslerian Poincare 
disc we cited above has a Yasuda-Shimada pedigree as well. Given that, it 
is desirable to have a more accessible and geometric proof of the theorem. 
The outcome of such an endeavor will be reported elsewhere. 

Finally, in §8, we give a rather preliminary discussion on the geodesies 
of our Randers metrics. Specifically, we invoke standard results about the 
projective Weyl and Douglas tensors to deduce that there is no coordinate 
system in which the geodesies appear as straight lines. Thus, as promised, 
our Randers spaces are not projectively flat. 



2. Randers spaces and the Yasuda-Shimada theorem 



Randers metrics were introduced by Randers in 1941 [Ra] in the context 
of general relativity. They play a prominent role in Ingarden's study of 



electron optics (see his treatment of the subject in [AIM]). Mathematically, 
in spite of the wide range of non-Ricmannian phenomena they are capable 
of producing, Randers spaces are Finsler spaces built from data that are 
quite familiar to all differential geometers: 

dx % ® dx J on a smooth n-dimcnsional 



• a Riemannian metric a 
manifold M, and 

• a differential 1-form b := bi dx 
drift 1-form. 



on M. This b is sometimes called a 
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Together these objects define a Finsler metric F in a simple way: 
F(x,y) := a(x,y) + 0(x,y), 

where 



<x(x,y) := J a y (x ) y V' 
/3(ic,y) := b l{x) y l . 



Here, a; stands for points on the manifold M, and j/ € T" X M denote tangent 
vectors based at x. Those tangent space coordinates y l typically come from 
the expansion y = y l -^~i in terms of a local coordinate basis. They can also 
arise, as y p , from the expansion y = y p e p in terms of an a-orthonormal 
frame field on M. Both scenarios are exemplified in §5. 

Due to the presence of (3, the Randers metric F := a + (3 is generally 
only positively homogeneous of degree one in y: F(x, cy) = cF(x, y) for all 
positive c. A Randers metric cannot be absolutely homogeneous [F(x, cy) = 
\c\F(x,y)] unless b = 0, in which case F is Riemannian. 

The fundamental tensor is formally analogous to the metric tensor in 
Riemannian geometry. It is defined as 



where we have used y 1 , y J as subscripts to signify partial differentiation. 
Almost by inspection, we have 



a 

Ei . — @i ~\~ hi. 

The fundamental tensor can then be expressed as 

Qij (j^ij ^j^j ■ 

Since f3(x, y) is linear in y, it cannot possibly have a fixed sign. The size 
of b therefore needs to be controlled if F is to be positive on TM \ 0. We 
also want the fundamental tensor to be positive definite. It turns out that 
both these properties hold if and only if 



|| 6 || := \lb, tt < 1, 
where 

V := cPlj. 

So, the drift 1-form b of Randers spaces must be required to have Riemann- 



ian norm strictly smaller than 1 everywhere. See |BCS or [ AIM | 
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The notion of flag curvature makes sense for all Finsler spaces. It is con- 
structed from the /i/i-curvature tensor of one's favorite Finsler connection. 
Happily, the resulting flag curvature is independent of which standard con- 
nection one is using, be it Berwald's, Cartan's, Chern's, or Hashiguchi's, 
just to name a few. Flag and spray curvatures will be reviewed in §6 rather 
than here, in order to effect a more streamlined exposition. 



The Yasuda-Shimada theorem [YS] says that for a Randers metric to 



have constant positive flag curvature K, the following four criteria are both 
necessary and sufficient: 

• The 1-form b is a Killing field of the Riemannian metric a. 

hi + hi = °- 

• The Riemannian norm of this Killing field must be constant, besides 
being globally less than 1. 

1 > II b II 2 := bi b % is constant. 



Second order covariant derivatives of b are to have the specific form 



' (a%k bj - a jk b^j 



h\j\k = K ( d ik b 

• The Riemann curvature tensor Rhijk of the metric a must be given 
by the special formula 

- K (l - || b |j 2 ) d hj d ik - K (d hj b t b k + d ik b h b^j + b h \j b~i\ k 

+ K (l - \\b\\ 2 ^ d h k dij + K {d h k h bj + dij b h bkj - b h \ k b^ 

+ 2 b h \i bj\ k . 

Here, the vertical slash (• • ')ju denotes covariant differentiation on M, taken 
with respect to the Levi-Civita (Christoffel) connection of the Riemannian 
metric d. 

In order to obtain a non-Riemannian Randers metric from these four 
criteria, the drift 1-form needs to be nowhere zero because it has constant 
length. Limiting our search to compact oriented manifolds M without 
boundary, we deduce from the Poincare-Hopf index theorem that the Eulcr 
characteristic x{M) must vanish. 

In two dimensions, the only candidate for M is therefore the torus. On 
this manifold, the most easily visualized Riemannian metric d is the one 
that gives the torus of revolution in M 3 . In that case, a straightforward 
calculation shows that every constant length Killing field b is identically 
zero. So there is no non-Riemannian Randers metric with constant positive 
flag curvature on the torus of revolution. 
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It is not obvious whether the above conclusion again results if we use 
other Ricmannian metrics a on the torus. If we are in the boundaryless 
category, there are telltale indications [^|]) that complete, positively 

homogeneous [as opposed to the more restrictive absolute homogeneity 
F(x,cy) = \c\F(x,y)] Finsler metrics of constant positive flag curvature 
can only be supported on manifolds homeomorphic to spheres. 

3. The Lie group S 3 and its Hopf fibration 

Our goal is to find a non-Riemannian Randers metric by 'solving' the 
four criteria stated in the Yasuda-Shimada theorem. This means we must 
produce a manifold M and, living globally on it, a special Riemannian 
metric a and a nonzero differential 1-form b that is appropriately coupled 
to a. As explained near the end of §2, the base manifold M must have 
zero Euler characteristic. Working in odd dimensions and staying in the 
boundaryless category automatically satisfies this topological constraint. 

The simplest compact boundaryless oriented 3-manifold is the Lie group 
S 3 , which is also a circle bundle over S 2 . We briefly discuss these properties 
for the sole purpose of setting some notation. 

• As a Riemannian manifold, S 3 is the standard unit sphere in Eu- 
clidean R 4 , whose points x have Cartesian coordinates (x°, x 1 , x 2 , x 3 ). 
In other words, S 3 can be characterized by 

\x\ := y/{x°) 2 + (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 = 1. 

• As a Lie group, S 3 consists of the unit quaternions 

x = x°l + xH + x 2 j + x 3 k 

in the space H of quaternions. Its group structure is isomorphic to 
SU(2). Multiplication between quaternions is non-commutative and 
is governed by the formal rules 

ii = 1, jj = 1, kk = 1, 
ij = +k, jk = +i, ki = +j, 
ji = ij, kj = jk, ik = ki, 

where the object 1 acts just like the real number 1. In view of the 
remarkable formula \xx\ — \x\ \x\, the unit quaternions are closed 
under the defined multiplication. Also, xx = \x\ 2 — xx, with 

x := x°l - x 1 ! - x 2 j - x 3 k. 

So the group inverse of any x £ S 3 is its conjugate: 
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As a circle bundle, S 3 has the celebrated Hopf fibration. For that pur- 
pose, it is best described as the set of 

(z, w) + ix 1 , x 2 + ix 3 ) 

in C x C satisfying \z\ 2 + \w\ 2 = 1. Incidentally, with this notation, the 
aforementioned isomorphism between S 3 and SU (2) reads 

(z, w) 

The unit complex numbers e lB S S* 1 act on S 3 by 

(z,w) i — ► (e w z, e t6 w) . 

The quotient manifold is the same as C 2 mod first the length and then 
the argument of an arbitrary complex number. So it is CP 1 which, by a 
suitable stereographic projection (see PJ), is precisely the complex manifold 
S 2 . Thus S 3 is a circle bundle over S . This bundle is non- trivial because 
S 3 is simply connected, whereas S 2 x S 1 is not. 

Using a curve of angles 9(t) with 6(0) — 0, one calculates the infinitesimal 
generator 

- (e w ^z,e^w) 
dt\t=o V / 

of the said 5* 1 action. The answer is 0'(0) times the following distinguished 

vector field 

E 1 := -x l d + x°di ~ x 3 d 2 + x 2 d 3l 

where abbreviates d x n . This E\ is therefore tangent to the S 1 fibres. It 
is complemented by 

E 2 := -x 2 d a + x 3 d 1 + x°d 2 - x 1 d z 

E 3 := -x 3 d Q - x 2 d 1 + x 1 d 2 + x°d 3 

to give a global orthonormal frame field on S 3 with 

[Ei, E 2 ] = — 2E 3} [E 2 , E 3 ] = — 2Ei, [E 3 , Ei\ — — 2E 2 . 

At any point x E S 3 , the values of Ei, E 2 , E 3 are respectively equal to the 
quaternion products ix, jx, and kx. This fact can be used to check that 
our frame field is right invariant. 

The natural dual of {Ei, E 2 , E 3 } is the following globally defined right 
invariant orthonormal coframe on S 3 : 



e 1 
e 2 
e 3 



-x 1 dx° + x°dx 1 - x 3 dx 2 + x 2 dx 3 , 
-x 2 dx° + x 3 dx 1 + x°dx 2 - x 1 dx 3 , 
dx° - x 2 dx 1 + x 1 dx 2 + x°dx 3 . 



-x 
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The standard metric on S 3 is thus 6 1 ® 9 1 + 9 2 <g) 9 2 + 9 3 ® 9 3 . Also, 
one has 



2 



Motivated by the treatment of the Hopf fibration in [GLP], we consider 
the following Riemannian metric on S 3 : 

a := e 2 9 1 <g> 9 1 + 9 2 <g> 9 2 + 9 3 <g> 9 3 , 

where e is a positive constant (typically different from 1). This a modifies 
the standard metric on S* 3 by introducing a dilation along the S 1 fibres. 
The Lie derivative of a can be calculated by using the Cartan formula 
£x6 = ixdO + dix® on 9 1 , 9 2 , and 9 3 . We find that every constant 
multiple of E\ is a Killing vector field of a, whereas E^, E3 are not Killing 
fields unless e happens to be 1. 

The Riemannian metric a admits its own orthonormal frame field 

ei := — Ei, e-i :— E2, &z '■= E3. 
e 

The natural dual consists of 

w 1 :=e0 1 , cj 2 :=9 2 , uj 3 := 9 3 , 

with 

2 2 
dui 1 = 2euj 2 A uj 3 , du) 2 = - uj 3 A U , dui 3 = - uj 1 A u? . 

e e 

Relative to this ct-orthonormal frame field, the Levi-Civita (Christoffcl) 

connection of a consists of a skew-symmetric 3x3 matrix of 1-forms u> q p . 

These are obtained by solving the structural equations duj p = ui q A u> q p and 

Lu pq = —ujqp. Here, ujqp means ui q s 5 sp , which is numerically the same as 

ujP. We find that 



/oj/ UJ 2 ijJ 3 \ ( -ELU 3 



UJn 1 LUn 2 




ywg 1 w 3 2 uj 3 3 J \-euj 2 (| - e) uj 1 
The curvature 2-forms and the Riemann curvature tensor are related by 
du q p ~uu q s Aw/ = i R* r „ uj t A oj s , 



where our convention for indices on R q p rs follows that in [ BCS |. Since we 
are in an orthonormal frame, R q p rs is numerically equal to R q p TS . Compu- 
tations give 

-R1212 = — e 2 , -R1213 = 0, -R1223 = 0, 

-R1312 = 0, -R1313 = — e 2 : -R1323 = 0, 

-R2312 = 0, -R2313 = 0, -R2323 = 3e 2 — 4. 
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All other components of the Riemann curvature of a are obtained from 
these by standard properties. Namely, R qpsr — —R qpr s, Rpqrs — —Rqprsi 
and the block symmetry R rsqp = R qprs - 



4. A privileged Killing field on S 3 

In §3, we discussed the Hopf fibration of 5* 3 and introduced the vector 
field Ei that is tangent to the S 1 fibres. The natural dual of this E\ is 
the 1-form 8 1 , also explicitly presented in §3. We pointed out that every 
constant multiple of E\ is a Killing vector (equivalently, every constant 
multiple of 1 is a Killing covector) field of the Ricmannian metric 

a := e 2 9 1 <g> 6 1 + 6 2 ® 6 2 + 9 3 <g> 6 3 

Motivated by this fact, let us stipulate the drift 1-form to be 

e 

where A is at the moment an arbitrary but nonzero constant. With respect 
to the metric a, this b has constant Riemannian length 




So the pair a and b do satisfy the first two criteria (among four) given by 
the Yasuda-Shimada theorem. Recall that those four criteria are necessary 
and sufficient for the Randers metric with data a, b to have constant positive 
flag curvature K. The purpose of this section is to use the remaining two 
criteria to determine the constants e and A in terms of K. 

We carry out our calculations in the moving frame {ei, 62, 63} and mov- 
ing coframe {u 1 , uj 2 , iv 3 } , which are o-orthonormal. The connection forms 
uo q p arc displayed in §3. We have 

h = - , b 2 = 0, b 3 = 0. 
e 

And the relevant covariant differentiation formulas are: 
b p \ g = {db p - b s uj p s ^ (e q ), 

bp\q\ r = (db v \ q - b s \ q uj p s - b p \ s uj q s ^ (e r ). 

Note that b p \ q is skew-symmetric on its indices because b is Killing. Con- 
sequently, the second covariant derivative b p \ q \ r is skew-symmetric in the 
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indices p and q. Straightforward calculations give: 

6i| 2 &l|3 N 
t>2\l b 2 \2 h\3 

\h\i h\2 h\3j 
and 

&1|2|1 = 0, 6i|3|i = 0, 6 2 |3|1 = 0, 

&1|2|2 = ~^ e j &1|3|2 = 0, 6 2 |3|2 = 0, 

&1|2|3 = 0) &1|3|3 = ~^ e > &2|3|3 = 0. 

Recall the third criterion given by Yasuda-Shimada. It reads 

bp\q\r = Tpq T , where Tpqr ■— K [dp r bq O-qr bp^j . 

We find that 

Ti2i - 0, T 131 = 0, T 231 = 0, 

^122 = — ^F, Ti 32 = 0, T 2 32 = 0, 

7 1 123 = 0, T133 = —^r, T233 = 0. 

Since the constant A is nonzero, the third Yasuda-Shimada criterion holds 
if and only if the dilation factor e is given by 

e = Vk. 

We have just determined that, for our purpose, the Ricmannian metric 
d on S 3 should be 

d := kb 1 ® e 1 + e 2 <g> e 2 + e 3 ® e 3 

This information updates the Riemann curvature tensor we calculated (at 
the end of §3) to 

#1212 = —K, #1213 = 0, #1223 = 0, 

#1312 = 0, #1313 = —K, #1323 = 0, 

#2312 = 0, #2313 = 0, #2323 = 3K — 4. 

Now examine the fourth (and last) criterion in the Yasuda-Shimada 
theorem. It requires the coupling between the drift 1-form b and the Ric- 
mannian metric d to be such that 

#(/prs Tqprs ; 
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where T qprs abbreviates the expression 

- K - || b \\ 2 ^j d qr a ps - K (a qr b p b s + a ps b q b r ^j + b q \ r b p \ s 

+ K (l - || b || 2 ) a qs a pr + K (a qs b p b r + a pr b q b s ^j - b q \ s b p \ r 

+ 2b q \ p b r \ s . 
Calculations give 

"71212 = —K, T~\21Z = 0, Tl22Z = 0, 

T\z\2 = 0, ^1313 = —K. 7i323 = 0, 

"72312 = 0, "72313 = 0, ^2323 = 4A 2 — K. 

Therefore the fourth criterion of Yasuda-Shimada is satisfied if and only if 

A = ±VK-l =: sVK-1 . 

In particular, the constant positive flag curvature K that one is striving for 
must be ^ 1. Anyway, the drift 1-form has now been determined. It is 

b := ±VK - 1 6 1 = ± 

It is somewhat amazing that the scaling multiple A (on b) asserts itself only 
at the -R2323 = 72323 stage. 

When K is equal to 1, the scaling multiple A vanishes and the dilation 
factor e is 1. The Randers metric in question then reduces to the standard 
Riemannian one that S 3 inherits from Euclidean R 4 . This case is not of 
interest to us because we want non-Riemannian Randers spaces. So let us 
impose the restriction 

K > 1. 



5. TWO EXPLICIT DESCRIPTIONS OF THE RESULTING RANDERS METRIC 

Let us recapitulate by giving explicit formulas for the Finsler functions 
of the Randers spaces we have just obtained. 

The first description is in terms of the non-holonomic frame {ei, e2, 63} 
and its coframe {u 1 , iv 2 , iv 3 }, both globally defined on the manifold M = S 3 . 
The frame is non-holonomic because each e p is a constant multiple of E p , 
and the latter have nonzero Lie brackets amongst themselves (see §3). 

Anyway, the Riemannian metric is 

a = lu 1 <g) uj 1 + uj 2 <g) lu 2 + lo 3 ® w 3 , 
and the drift 1-form is 
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Expanding arbitrary tangent vectors as 

V = y 1 ei + y 2 e 2 + y 3 e 3 , 
we see that the Finsler function is given by 



F(x, y ) = Viy 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 ± 



K - 1 




K 



It is remarkable that there is only an implicit dependence on the position 
x. We will revisit this deceptively simple formula in §7. 

The second description of our Randers metric is in terms of natural 
coordinates. For this purpose, we used the following parametrization of S 3 : 



where c has the value +1 when dealing with the right hemisphere, and the 
value —1 when dealing with the left hemisphere. Obviously, the y here 
no longer denotes generic tangent vectors on S 3 . Instead, it is part of 
the collection x, y, z that one typically uses when discussing Cartesian 
coordinates in Euclidean R 3 . 

This parametrization has the advantage that it imposes no restriction 
on x, y, z, and can be visualized as follows. Consider for example the right 
hemisphere of S* 3 , centered at the origin O := (0,0,0,0) of M 4 . Place a 
hyperplane R 3 tangent to the sphere at the East Pole (1,0, 0, 0). Points on 
this tangent hyperplane are of the form (l,x,y, z). Now we assign coordi- 
nates to any given position P := (x°, x 1 , x 2 , x 3 ) on the right hemisphere of 
S 3 . Draw the straight line segment from the origin O to P and prolongate 
that until it intersects the tangent hyperplane. The point of intersection, 
being on that hyperplane, will have the form (l,x,y, z). Just declare the 
natural coordinates of our point P to be (x,y,z). Note that points along 
the equator of the right hemisphere correspond to points at infinity on the 
tangent hyperplane. A similar story holds for the left hemisphere. 

The Riemannian metric a and the drift 1-form b are 



(^X ^ X ^ X j X ^ — 



1 



= (c,x,y,z) 



y/1 + x 2 + y 2 + z- 



a= Ke 1 ®e 1 + o 2 ®& 2 + e 3 ®e 3 



b = ±vk - i e 1 . 



Simple computations give: 
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+ cdx — zdy + ydz 



e 2 




e 3 



l + x 2 +y 2 + z 2 ' 
— ydx + xdy + cdz 



1 + x 2 + y 2 + z 2 
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Denote generic tangent vectors on S 3 as 

d d d 
ox Oy oz 
Then the Finsler function for our Randers space is 

F(x, y, z; u, v, w) — a(x, y, z; u, v, w) + (3(x, y, z; u, v, w), 

with 

\l K(cu — zv + yw) 2 + (zu + cv — xw) 2 + {—yu + xv + cw) 2 

CL = — 

1 + x 2 + y 2 + z 2 



±V K — 1 (cu — zv + yw) 
1 + x 2 + y 2 + z 2 

As a reminder, c = +1 for the right hemisphere, and c = — 1 for the left 
hemisphere. Some Maple programming will be carried out on this formula 
in the Appendix. 

The Riemannian norm of the drift 1-form is 



It is given by this value whether one is using the a-orthonormal frame or 
natural coordinates. Since K ^ 1 (note: K = 1 corresponds to the standard 
sphere S* 3 ), we see th a t \\b\ \ is strictly less than 1. This is both necessary 



and sufficient ([ AIM |, [BCSQ for F to be positive and strongly convex on 



the entire slit tangent bundle TM \ 0. 



6. Spray curvatures and flag curvatures 

The flag curvature is defined much like the sectional curvature of Rie- 
mannian geometry. One begins with a connection, differentiates appropri- 
ately, and then performs some standard contractions. Let us first give a 
synopsis of the procedure. 

The connection in use (from a family that contains the handful of name 
brand ones) gives rise to a /i/i-curvature tensor J? * which in the Riemann- 
ian case is precisely the familiar Riemann curvature. The precise formula of 

kl in term s of, s ay, the Chern connection, has been given a pedagogical 



treatment in | BC£ |. Since that specific formula does not concern us in this 
paper, we shall omit it. 

At any point i on M, one associates with each flagpole — a nonzero vector 
y in T X M — the spray curvature 

k\ : ,rii; u ,/. 
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For those familiar with [ BCS| , this K l k is F 2 times something called the 



'predecessor of the flag curvature.' It is actually the first among two spray 



curvatures (|B3|, M). The second one, using the version described in 



and multiplying by F, is equal to the /iu-curvature h Pj % u of the Berwald 
connection, but is irrelevant to the purpose of this paper. 



Our notation K l k is the same as Berwald's jBl ], Rund [H uses H l k 
instead. The spray curvature is robust enough that it does not depend on 
our choice of connection. However, the above definition is often not an 
efficient way to compute the spray curvature because one has to obtain the 
tensor R/ kl first. 

A word of caution about notation. For us, Kn~ shall simply mean 



gijK 3 k , It can be shown (see BCS]) that Ki k is symmetric in 
i, k, and y l Ki k , Kn~y k both vanish. Rund did not subscribe to 
this system. He defined his Hn~ by first tracing the indices j 
and I on ^[{H° l ) y k — [H 3 k )yi], forming an intermediate H kl and 
then taking the y-partial (H k ) y i ■ See pages 129 and then 125 of 
(rJ. As a result, his Hi k was not symmetric and y l Hi k did not 
vanish for him. 

The flagpole y, together with any vector V <G T X M transversal to y, 
specifies a flag based at the point x € M. The flag curvature of the resulting 
flag is the quantity 

V 1 K j V k 

K(y,V):= 



9(y,y)g(v,v)-[g( y ,v)f 



where g := g^dx 1 dx^ is the fundamental tensor defined in §2. The 
Finsler metric is said to have constant flag curvature if K (y, V) has the 
same constant value K for any choice of y and V . 

Specialize now to the case of constant flag curvature K. Euler's theorem 
for homogeneous functions implies that g(y, y) — F 2 (x, y). Let's substitute 
this into the above formula and rearrange it into the form 

V* K ik V k = K F 2 {g(V, V) - [g(£, V)] 2 }, 

where I means y/F. As such, this statement makes formal sense even if V 
is not transversal to y. Since both sides are given by symmetric bilinear 
forms, a standard polarization identity gives 

U l K ik W k — K F 2 {g{U, W) - g(£, U)g{£ 7 W)}. 

In particular, 

Kik = K F 2 (g ik - F yi F yk ) . 

Here, we have used the identity g S j£ s = F y j , again a consequence of Euler's 
theorem. Raising the index i with the inverse matrix g lJ of the fundamental 
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tensor, we get 

We shall use this criterion as the characterization of constant flag curvature. 

Next, we describe a computationally friendlier way to access the spray 
curvature K l k . The history behind this better approach goes all the way 



back to Berwald |B1| in his study of path spaces and projective geometry. 
See Rund Q for a careful exposition, and references therein. 
One begins with the geodesic spray coefficients 

where 

i t i 1 ( dgij dg jk dg k[ 



' jk ' y 2 V dx k dx l dxi 
are the fundamental tensor's formal Christoffel symbols of the second kind. 
Our definition of G l here is | times that given in the book [ BC£[|. As a 



result, the present definition agrees exactly with the one used in [AIM] and 
Chapter IV formula (6.3) in tells us that 



K\ = 2{G% k - y>{G l ) xlyk - {G%{G%, + 2&{G i ) 



y J v 



For ease of exposition, let us refer to this as Berwald's formula. It is an 
instructive exercise to manipulate the original expression for the spray cur- 
vature, namely yi iZ * kl y l , into the one we just displayed. See §12. 5B in 



[BCS| if guidance is needed. 

The above formula expresses the spray curvature in terms of partial 
derivatives of the geodesic spray coefficients. As such, it is quite useful 
for machine computations of the spray curvature. The Maple codes that 
implement this sort of computation are given in the Appendix. 

The said formula is also invaluable for calculating spray curvatures by 
hand. However, its efficient implementation for this purpose often requires 
some additional maneuvers as described in |^l). To get us set up for §7, we 
specialize this technique to the Randers setting and present it here. 

The geodesic spray coefficients of every Randers metric is equal to that of 
the underlying Riemannian metric, plus a perturbation term. Symbolically, 
we write 

G l = G l + C- 

Here, 
where 



7 jk ■= a 



_ a 1 / daij _ dhjk dau 
2\dx k dx l dxi 
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are the Christoffel symbols of the second kind of the Riemannian metric 
a. The specific form of the perturbation term £ ! has been worked out in 
[ BCS |, but shall not concern us until §7. 

Substitute the above G z into Berwald's formula for the spray curvature. 
After covariantizing (to be explained below) all x-partial derivatives of 
and using Euler's theorem for homogeneous functions whenever appropri- 
ate, we get 



This formula will be our centerpiece in §7. The first term on the right-hand 
side, K l kl is the spray curvature of the Riemannian metric a. It can readily 
be calculated in the S 3 examples because the Riemann curvature Rjiki of 
a is already at our disposal (§4). 

The remaining group of terms on the right-hand side involves horizontal 
covariant differentiation on the slit tangent bundle TM \ 0. Let us explain 
how those are carried out. The geodesic spray coefficients G l are first used 
to define 

and then 

5 d ^ d 
The vector fields -J-r are declared to be horizontal the horizontal lift of J~ 

OX J _ ox° 

to TM \ 0. And the iV l - are said to have produced a nonlinear Ehresmann 
connection on the slit tangent bundle. In order to achieve the sought hor- 
izontal covariant differentiation on we horizontally lift the Levi-Civita 
(Christoffel) connection of a, and then let the resulting object act on Q . 
This may seem abstract but is operationally quite simple: 

Incidentally, note that 

This complements nicely with the fact that the nonlinear connection is 
given by the first y-partial derivatives of G l . Thus 

{C)xi = C\j + (C) y '(G l )yi - C\G l ) y i y] . 

The use of this formula is what we mean above by covariantizing the im- 
partial derivatives of Q . We did that to make each term inside the expres- 
sion {• • • } manifestly tensorial. Of course, the quantity {• • • } as a whole 
is already tensorial, whether we covariantize or not. There is no need to 
covariantizc the y-partial derivatives because they already transform ten- 
sorially (see BCS| for an exposition). 
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7. Direct proof of constant positive flag curvature 

As we mentioned in the Introduction, the two published proofs (| 
and Q ) of the Yasuda-Shimada theorem both contain arguments that we 
find difficult to follow. Given that, it is prudent to prove directly that our 
Randers metrics have constant positive flag curvature K. The numerical ev- 
idence documented in the Appendix helps sustain such a bare-hands' proof. 
Strategically, we compute both sides of the constant curvature criterion 



K p „ = KF 2 [S p - KtF, 



F 



and ascertain that they are indeed equal. The calculations rely on Berwald's 
formula, and will be done with the help of an a-orthonormal frame instead 
of natural coordinates. The resulting proo f is geometrical and is distinctly 
different from the approaches used in [YS], jMj. It presents a new perspec- 



tive on the subject, and is worthy of further development. 



7.1. The right-hand side of the constant curvature criterion. 

Let us begin with the formula for F in terms of the a-orthonormal frame 
{ei,e2,e3} on S 3 . This was first given in §5. For ease of presentation, we 
relabel the tangent space coordinates y 1 , y 2 , y 3 as u, v, w, respectively. 
These coordinates arise from the expansion of arbitrary tangent vectors y 
in terms of the basis {ei,e2,e3}. The cosmetically altered formula for F 
now reads 



F(x, y) — \/u 2 + v 2 + w 2 + 5\j ^ ^ ^ u =: a + (3 . 

Here, s — ±1 keeps track of the two choices of the drift 1-form. See the 
end of §4. These choices of sign have nothing to do with the right or left 
hemisphere of S 3 . The latter are associated with the parameter c of §5. 
Almost by inspection, we have 



F, 



a V K 



K - 1 



a 

F =- 

a 

Using these, the right-hand side 



KtK :=KF 2 ( S^-'Cf, 



F 
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of our constant curvature criterion is readily computed. We find that 



Kt\ Kt 2 



Kt 2 



K t% K r 3 



is equal to 



V 



£ KN 2 + w 2 

a v 

^Ku\j-FK\js 

a 

— Kuw — F Kws 



K-l 
K 

K-l 
K 



a 

F K(F — ^ 
— — Kwj 



a 

— Kv\n 



FK(F~^)J 



7.2. Perturbation terms in the geodesic spray coefficients. 

The left-hand side of our constant curvature criterion is the spray cur- 
vature K p r . These nine components can be calculated using Berwald's 
formula. To this end, we start with the geodesic spray coefficients G\ As 
we mentioned in §6, the ones used in this paper, unlike those of [BCS], 
already have the factor of ^ built in. It is shown in the same reference that 
for Randers metrics, 

where G l := \ 7* - fc y^y k are the geodesic spray coefficients of the underlying 
Riemannian metric a. 

The perturbation terms £ l transform like the components of a tensor. 
They do so because 

showing that ( arises from the difference of two connections. More ab- 
stractly, consider the vector bundle TM that sits over M. Using the pro- 
jection map 7T : TM \ — > M, we can pull that back to obtain a vector 
bundle ir*TM that sits over the manifold TM \ 0. This just means that 
over each point (x, y) € TM \ 0, we have erected a copy of T X M . Our Q 
then transforms like a section of this pulled-back vector bundle. 

The coordinate bases {d x i}, the Riemannian metric a, together with the 
field of a-orthonormal frames {e p }, can all be transplanted to the fibres of 
ir*TM (note: not to TM \ 0). As on M, the transplants of {d^} and {e p } 
are related through the 3-bein w z and its matrix inverse 



' d,. 



dx\ 



Each and is a function on M. In other words, they depend on x 
only. 

The tensorial property mentioned above allows us to transform the co- 
ordinate description Q into the orthonormal frame description 
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With this in mind, and the factor of | emphasized earlier, an expression 
derived in BCSfl is transcribed. It gives 

Jp> y«y r 



(a pq y r 



~a pr y q ) 



£P (yl V _ f b qj 



Here, £ p := y p /F is to be distinguished from l p :— y p /a. We shall show 
that these ( p components are actually quite simple. 

In §4, we determined the intermediate constants e and A. These can be 
used to update the Levi-Civita (Christoffel) connection io p of the Riemann- 
ian metric a, the components of the 1-form 6, and its covariant derivatives 



A quick glance at §3 and §4 tells us that the updated versions read: 



"l 1 


-I 2 






uj 2 2 


UJ 2 3 


Wg 1 


iv 3 2 


UJ 3 3 





Kuj 3 



Kco 3 




K " 2 {ve-^)" 1 



'Ku? 




together with 



and 



An 

\h\i 




0, b 3 



&2|2 b 2 \3 
h\2 h\3 









1 








Since the present b q \ r is skew-symmetric and most of its components are 
zero, the perturbation term £ p simplifies drastically to: 



a 5 VK~^T {a p3 y 2 - ~a p2 y 3 + £ p (y 3 b 2 - y 2 6 3 )} 



Furthermore, we are in an a-orthonormal frame, so b r := a rs b s = 8 rs b s is 
numerically the same as b r . Given that 62 and 63 are both zero, the terms 
multiplying l p drop out, and 



In other words, 



( p = asVW^T{S p3 y 2 - S p2 y 3 } 
0, 



c 1 

( 2 = -sV K - 1 wa, 
C 3 = +WK - 1 va. 
Recall that in the a-orthonormal frame description, the quantity a is simply 



w 2 . So the components £ p contain no explicit dependence on 
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x. However, u, v, w are y p , which can be expressed as v p i y l . Thus the 
dependence on x is there, albeit implicitly through v p i . 

7.3. The spray curvature of the Riemannian metric. 

In §6, we described a refined understanding of Berwald's formula. It is 
a case study of the general strategy discussed in ||Sl|] . The essence is that 
given the decomposition G l — G l + C , the spray curvature splits as well: 



Each side of this equation transforms like a section of the pulled-back tensor 
bundle ir*TM <g> ir*T*M, which sits over TM \ 0. 

Contract both sides of the above equation with v p i u r k . This converts its 
expression from the coordinate to the a-orthonormal perspective: 

K p r = K> r + {2C p lr ~ y 9 (C p k V - (C P MCV + 2C 9 (C p ) y ^} • 

The first term on the right-hand side is the spray curvature of the Rie- 
mannian metric a. Since the Riemann curvature tensor Rq prs of a has been 
determined in §4, we can use that information to compute the said spray 
curvature: 

Keep in mind that we are in an a-orthonormal frame, so the numerical 
values of tensor components are unaffected by raising or lowering any index, 
as long as it is done using a. Fairly routine calculations tell us that the 
nine spray curvatures 

' K\ K\ K\ 

k\ k\ k\ 
k\ k\ k\ 

are equal to 

'if(v 2 +w 2 ) -Kuv -Ku\n 

-Kum Ku 2 + (4 - 3K)w 2 (3K - 4)vw 
-Kuw (3K - 4)vw Ku 2 + (4 - 3K)v 2 

The symmetry here is expected because K p r := a ps K sr has the same nu- 
merical value as K pr , and the latter is symmetric by general principles. In 
contradistinction to this, the indices on K p r are raised and lowered with 
g rather than a. Since the basis we are using on the fibres of tt*TM is 
not g-orthonormal, we do not expect K p r to be numerically equal to the 
symmetric K pr . 
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7.4. A simplified formula for horizontal covariant derivatives. 

In the refined version of Bcrwald's formula, we have the terms 

These involve both the y-partial derivatives of £ p as well as its horizontal 
covariant derivatives. The purpose of this subsection is to derive a simplified 
formula for these horizontal covariant derivatives, one that will facilitate 
their computation in §7.5. 

In §6, we explained how the horizontal covariant derivatives are to be 
carried out in the context of natural coordinates. Namely, 

/■i /•% _i_ A ~i 

|fc ~ So 1 lk ' 

Note that 7^ ■ dx^ are the connection forms with respect to the coordinate 
basis. When these connection forms are pulled back to TM \ and evalu- 
ated on the horizontal lifts (of d x k to TM \ 0) 

5 a ~ H _d_ 

8x k dx k k dy % ' 
they see only the d x k part but not the d y i part. The above expression for 
can then be recast as 

C lk = {dC + C i (f y ^)}(^) • 

In the current subsection, we are using the basis {e p ®u) r } (instead of the 
coordinate one) on the pulled-back bundle n*TM®ir*T*M. The analogous 
formula for ( p , is 

( P lr = {d( p + C s ^/}(e,) , 
where tu s p are the connection forms displayed in §7.2. The vector field e r 
is the horizontal lift of e r , and has the formula 

e r := e r — y s lj s p (e r ) d yP . 

This is not unreasonable, given that in the horizontal lift of d x k , the non- 
linear connection can be re-expressed as follows: 

N\ = i\ k y l = y l (7*y (a,*). 

Note, however, that we have yet to define e r , and for that matter d x k, 
as objects on TM \ 0. Here are the definitions: 

• The 1-forms dx k are pulled-back to TM\0 and given the same name. 
They are then complemented by the dy k (with y k coming from y k d x k) 
to form a coordinate basis for the cotangent bundle of TM \ 0. The 
natural dual of this basis is denoted {d x k , d y k}. The objects in this 
basis are local vector fields on TM \ 0. They are the ones that enter 
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the definition of . It is simply an abuse of notation to employ the 
same symbol d x k for the vector fields here and those on M. 

• In the same spirit, we pull back the I-forms uo r to TM \ and retain 
the same name. These are complemented by the dy r (with y r arising 
from y r e r ) to give a non-holonomic basis. Denote the natural dual 
of this basis as {e r ,d yr }. Just like the d x k case, we are abusing the 
notation when we call these e r and those on M by the same name. 
Nonetheless, we now know what e r means on TM\0, and the formula 
for the horizontal lift e r makes good sense. 

Do not confuse the e r and d x k we just defined with the transplants (see 
§7.2) that live on the pulled-back vector bundle tt*TM, even though they 
share the same notation merely for the sake of economy. 

One can't help but wonder how e r is related to d x k as vector fields on 
TM \ 0. To find out, take an arbitrary differentiable function / on TM \ 
and let e r act on it. By the chain rule, 

(df)(e r ) = {(d x kf)dx k + (d yk f)dy k } {er). 

• Note that y k = y s u s k gives dy k = u k dy s + y s du k . It is to minimize 
confusion in this statement that we have chosen the font y s throughout 
the paper, when dealing with orthonormal expansions of y. 

• The relation dx k — u k uj q is valid on TM \ because it holds on M 
and the pull-back operation preserves algebraic statements. 

• Since {e p , d yP } is by definition the natural dual of {iv p , dy p }, one must 
have uj g (e r ) = S q r and dy s (e r ) = 0. 

Substituting these observations into our chain rule statement, and removing 
the test function / afterwards, we obtain the somewhat surprising formula 

e r = u k d x k + y s du k (e r ) d y k on TM \ 0. 

In contradistinction to that, a similar calculation gives 

dyr = u k d y k on TM \ 0. 

Here, one does have to invoke (du s k )(d y r) = 0. This holds because u k is a 
function on M, so its differential is a linear combination of dx^ and hence 
of uj q . Such a relation is preserved under pull-back, and Lo q {d y r) = 0. 
Let us return to the horizontal covariant derivative 

e ir = {dC p + C"/}(e r ) , 

with e r := e r — y s w s *(e r ) d y t. This can stand two more reductions before 
we put it to use in §7.5. 

• On M, the connection forms ut s p are linear combinations of the ou q , 
which arc in turn linear combinations of the dx l . This remains so 
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under the pull-back to TM \ 0, so u> s p (d y h) — 0. Since d y t — u t k d y k, 
we must have <jj p (d y t) = as well. Therefore 

W /(e r )=w/(e r ). 

• At the end of §7.2, we deduced that the components ( p depend 
only on u, v, w, namely y q . Formally taking the differential yields 
d( p = (( p ) y <idy q . On the other hand, linear algebra says that d( p = 
(d( p )(e q ) uj q + (d( p )(d yq ) dy q . Comparing the two statements gives 
(d( p )(e q ) = 0, which in turn allows us to conclude that 

K p )(e,) = -y s ^( er )(C P V. 

The two reductions above lead to the following: 

CV = {C S "/-(C P )y' y s LoJ}(e r ). 

We shall use this formula in §7.5 to compute horizontal covariant deriva- 
tives. By the way, in the context of [BCS], the lo p here happens to be the 



Chern connection forms of the Riemannian metric a. 

7.5. Derivatives of the perturbation terms. 

In this subsection, we tabulate the requisite y r -partial derivatives of £ p . 
Then we use the concluding formula of §7.4, together with lo p from §7.2 
and u> s {e r ) = 8 s r , to calculate the horizontal covariant derivatives that we 
need. Recall that 

C 1 = 0, C 2 = - V K - 1 wa, C 3 = 1 vd, 

with 

a := \J u 2 + v 2 + w 2 . 

Here, u = y 1 , v = y 2 , w = y 3 and, from §7.1, s = ±1. Thus all y r -partials 
of are zero. On the other hand, £ < r do not necessarily vanish because 
C 2 and C 3 are involved in the calculations. 

Let us list the results. We introduce some abbreviations for ubiquitous 
quantities in order to reduce clutter. 



With := — y/K-1 



(C 2 ) y i =-0(uw) (C 3 V =+0(uv) 



a 

i 

(C 2 )y^ = (VW) (C 3 )y2 = +0 (U 2 + 2V 2 + W 2 ) 

(C 2 )y^ = (u 2 + v 2 + 2w 2 ) (C 3 ) y3 = +0 (vw) 



Jy 1 

'2\ 
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(C 2 )yly2 = +* (UVW) (C 3 W = +* "(U 2 + W 2 ) 

(C 2 ) yly3 = -* u(u 2 + v 2 ) (C 3 ) yly3 = * (uvw) 

(C 2 ) y2y2 = -* w(u 2 + w 2 ) (C 3 ) y2y2 = +* v(3u 2 + 2v 2 + 3w 2 ) 

(C 2 ) yV = -* v(u 2 + v 2 ) (C 3 ) y2y3 = +★ w(u 2 + w 2 ) 

(C 2 ) yV - -* w(3u 2 + 3v 2 + 2w 2 ) (C 3 ) yV = +★ v(u 2 + v 2 ) 



With 9 := a s VK-l \[K : 

CV = o C 2 n = o c 3 n = o 

C 1 | 3 = -^W C 2 | 3 = C% = +^U 



Lastly, 



With * := - s y/K-1 \F~K : 



(C 2 |2 ) i = +* (2u 2 + v 2 + w 2 ) (C 3 |3 ) i = +* (2u 2 + v 2 + w 2 ) 

(c 2 |2 ) y2 =+*M (c 3 |3 )" 2 =+*M 

(C 2 |2 ) y3 =+*(uw) (C 3 |3 ) y3 = + *( U W) 
These two columns are equal because (C 2 2 ) = (C 3 ^)- Also, we have 

(c i |a ) yi = -*(uv ) (cg yi = -*(uw) 

(C 1 |2 ) y2 =-*(u 2 + 2v 2 + w 2 ) (^3)^=-*^) 
(C 1 |9 ) ys =-*(vw) (C 1 |3 ) y ^-*(u 2 + v 2 + 2w 2 ) 

7.6. The spray curvature of our Randers metric. 

According to §7.3, the spray curvature of the Randers metric 

F := a + (3 

has the structure 
where 

£? r := 2^ r - y"(C P lg ) y r - (C P MC"V + 2C 9 K P W- 

The term K p r represents the nine spray curvatures of the underlying Ric- 
mannian metric a. We have already calculated those in §7.3. 
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Using the derivatives listed in §7.5, the above quantities £ p r are com- 
puted. We find that: 



e\ = 


+*u(v 2 + 


w 2 ), 










E\ = 


-*v(u 2 ), 












e\ = 


-*w(u 2 ); 












e\ = 


-*v(2u 2 - 


f v 2 + 


w 2 )- 


(K- 


l)uv, 




f 2 - 


+*u(2u 2 - 


f v 2 4 


- 2w 2 ) - 


f (K 


-l)(u 2 - 


h4w 2 ), 


£\ = 


— ft uvw — 


A(K- 


- l)vw; 










-*w(2u 2 


+ v 2 H 


-w 2 ) - 


(K- 


- l)uw, 




F 3 - 


— W|b uvw — 


A{K - 


- l)vw, 








F 3 - 

C 3 — 


+*u(2u 2 - 


f 2v 2 


+ w 2 )- 


f {K 


-l)(u 2 - 


h4v 2 ). 



Adding these £ p r to the corresponding X p r of §7.3, we obtain the spray 
curvatures K p r . After that, it is a matter of routine algebra to make sure 
each K p r is equal to the Kr p r calculated in §7.1. Such is indeed the case. 
This means that the Randers metric in question satisfies the criterion 

K p r = KF 2 (s p r -^F yr y 

Therefore it has constant flag curvature K and our proof is complete. 



8. Discussion 

Given any Finsler manifold (M,F), its geodesies are paths in M, just 
like the Riemannian case. The defining equation for geodesies with constant 
Finslerian speed is 

x i +i i jk x j x k = 0. 

In terms of the coefficients G l of §6, this system of second order quasi-linear 
equations can be rewritten as 

x l + 2 G l = 0. 

This is one reason the G l are called the geodesic spray coefficients. 

A Finsler manifold is said to be projectively flat if M can be covered 
by privileged coordinate charts in which the above quasi-linear system of 
ordinary differential equations becomes a linear system. Projective flatness 
has been characterized by Douglas Jd| in his study of path spaces. When 
specialized to the Finsler setting, his result states that: 
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A Finsler manifold (M, F) of dimension ^ 3 is projectively flat if 
and only if its projective Weyl and Douglas tensors both vanish. 
A slightly modified characterization, due to Berwald [|B2| , applies to di- 
mension 2. There, the vanishing of the Weyl tensor is replaced by another 
criterion. See page 144 of Rund |f|], or for a detailed account and 
further references. 

The projective Weyl tensor is to be distinguished from the conformal 
Weyl tensor in Riemannian geometry. The latter is not defined in dimension 
2 and always vanishes in dimension 3. It turns out that the vanishing of 
the projective Weyl tensor (which has four indices) is equivalent to the 
vanishing of its reduced version W\. This is explained on pages 141 and 
142 of §§. 

The tensor W\ is defined in terms of the spray curvature and its deriva- 
tives. To that end, one first constructs the scalar quantity 

n — 1 



Comparing with §7.6 of [BCS|, we see that this scalar A can be interpreted 
as F 2 times the average of n— 1 appropriately chosen flag curvatures. Using 
the reduced projective Weyl tensor is defined as 



W\:=K\-&5\-- 
n 



A result of Matsumoto's (see [AIM]) says that the vanishing of W\ is 
both necessary and sufficient for the flag curvature to depend at most on 
the position x and the flagpole y S T X M, but not on the transverse edges 
(denoted V in §6). In particular, since our Randers metrics on S 3 have 
constant flag curvature K, we deduce that W\ must vanish for all these 
examples. This can also be verified by a direct computation. 

The Douglas tensor D l - kl (denoted B % - kl in Q) is defined as follows: 

D l jkl ^ {G^yjykyl 

— {5 l j(G h )yhykyl 4~ <5 2 fc ( G ^ ) y fl J, ! J, j "f <T J ( G ^ ) y h y j y k ^ 



n 

' y l (G h ) yh 



1 



y"yj y^y L 



In | BM |, Basco and Matsumoto proved that a Randers space has vanishing 
Douglas tensor if and only if the drift 1-form b is closed. For our examples, 
b is a constant multiple of 6 1 , see §4. Since dQ 1 = 29 2 A 9 3 / 0, we 
conclude that the Douglas tensor does not vanish for our Randers metrics. 

Now we know that each member among our family of Randers metrics 
on S 3 has vanishing projective Weyl tensor but nonzero Douglas tensor. 
Therefore, by Douglas' theorem, none of them can be projectively flat. 
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Appendix A. Some numerical evidence 

The purpose of this Appendix is to apply Maple to the Randers metric 

F(x, y, z; u, v, w) = a(x, y, z; u, v, w) + f3(x, y, z; u, v, w) 

in natural coordinates, with 

y/ K(cu — zv + yw) 2 + (zu + cv — xw) 2 + (—yu + xv + cw) 2 
a= 1 + x 2 + y 2 + z 2 ' 

±V K — 1 (cu — zv + yw) 
f = 1 + x 2 + y 2 + z 2 ' 

Recall that c = +1 for the right hemisphere of S 3 , and c = —1 for the left 
hemisphere. 

In Finsler geometry, it is not uncommon for simple formulae to quickly 
mushroom into unmanageable expressions. We find that machine compu- 
tations have consistently extracted useful information and insights, to the 
point that meaningful followup questions can be asked. We value every op- 
portunity to cultivate this synergy between Finsler geometry and modern 
computing. It is hoped that by producing the Maple codes here, we are ini- 
tiating other geometers into a fruitful aspect of experimental mathematics. 

Here is our plan: 

• We first use Maple to calculate the spray curvature a la Berwald: 

K\ = 2{G l ) x k - y j (G l ) x3y k - {G l ) y] {G\ k +2&{G l ) y]yk . 

• After that, we ask Maple to check whether F satisfies the characteri- 
zation of having constant flag curvature K: 

K\ = KF> ^5\- V -F y ^. 

A.l. The Finsler function in natural coordinates. 

> P:=c*u+y*w-z*v; 

> Q:=c*v+z*u-x*w; 

> R:=c*w+x*v-y*u; 

Let us restrict our attention to the right hemisphere of S* 3 . So 

> c:=+l; 

Define, for lack of a better name: 

> dcn:=l+x"2+y~2+z"2; 
Then 

> alpha:=(l/dcn)*sqrt(K*P"2+qT2+R"2); 

> beta:=(l/den)*sqrt(K-l)*P; 



2,x 
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Here, we have simply chosen the + sign in the drift term (3. A moment's 
thought shows that there is no loss of generality in doing so. The Finsler 
function F and its associated Lagrangian L are: 

> F:=alpha+beta; 

> L:=(1/2)*F~2; 

Use of the semi-colon instructs Maple to display the input formulas. 

A. 2. The covariant form of the geodesic spray coefficients. 

By that, we are referring to the quantities 

Gi -.= \ 7yfe y J y k = \ (.</,.,.,•• - 9jk,x* + 9ki,xi) v 3 y k - 

Raising the index with the inverse of the fundamental tensor gives the 
contravariant form G l . That will be carried out in §A.4. In terms of the 
Lagrangian L, it is not difficult to show that 

= 2 (Lyixi y° ~ L x i ) . 

Our natural coordinates x 1 , x 2 , x 3 are denoted x, y, z, and the induced 
tangent space coordinates y\ y 2 , y 3 are denoted u, v, w. The covariant 
form of the geodesic spray coefficients are named (by us) gl, g2, <?3 in 
Maple. Thus 

> Lx:=diff(L,x): 

> Ly:=diff(L,y): 

> Lz:=diff(L,z): 

> gl:=(l/2)*(u*diff(Lx,u)+v*diff(Ly,u)+w*diff(Lz,u)-Lx): 

> g2:=(l /2)*(u*diff(Lx,v)+v*diff(Ly,v)+w*diff(Lz,v)-Ly): 

> g3:=(l /2)*(u*diff(Lx,w)+v*diff(Ly,w)+w*diff(Lz,w)-Lz): 
Use of the colon instructs Maple to suppress the computed answers. 



A. 3. The inverse of the fundamental tensor. 

The inverse g lJ of the fundamental tensor gij (§2) is needed in order to 
obtain the contravariant form G % of the geodesic spray coefficients. This 
inverse has a standard formula 



F a ' 



F 2 



F 2 { 



t V + £ J b' 



where t 



a lJ bj , and a lJ is the inverse of the Ricmannian metric 



; . For a pedagogical derivation, one can consult [BCS 



In our natural coordinates, the matrix of the Riemannian metric a is 

' K + z 2 + y 2 —Kcz + cz — xy Key — xz — cy^ 
— xy Kz 2 + 1 + x 2 —Kyz 



1 

den 2 



-Kcz 



i Key — xz — cy 



+ 1 + 
-Kyz 



Ky 2 



1 
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with den :— l+x 2 +y 2 + z 2 . The inverse of dij has been computed elsewhere 
using Maple. We assign it the name Aij in our codes: 

> All:=(l/K)*((x~2+l)*(x~2+K*z~2+K*y~2+l)): 

> A12:=(l/K)*(x~3*y+xV3*K+x*y-c*z*x~2-c*z 

+K*c*z*x~2+K*c*z+K*y 5 V2*x): 

> A13:=(l/K)*(K*y~2*z*x+x~3*z+x*z~3*K+x*z 

-K*c*y*x'2-K*c*y+c*y*x~2+c*y) : 

> A21:=A12: 

> A22:=(l/K)*(y2*x"2+K*x"2+(K-l)*2*x*y*c*z 

+y"4*K+z"2+K+y~2*K*z"2+2*K*y"2): 

> A23:=(l/K)*(x~2*z*y+(K-l)*c*z~2*x-(K-l)*x*y~2*c 

+K*y'3*z+K*y*z'3-z*y+2*y*K*z) : 
A13: 
A23: 



> A31 

> A32 

> A33 



<l/K)*(z~2*x~2+K*x~2-(K-l)*2*x*y*c*z 



+y"2*K*z~2+z"4*K+2*K*z~2+K+y"2): 

As we have previously decided, c = +1 because we want to focus on the 
right hemisphere. The (covariant) components 6j of our drift 1-form are 



l + x 2 + y 2 + z 2 



(c, -z, +y). 



Denote its contravariant components b % := a 1J bj in Maple as Bl, B2, B3. 
On the right hemisphere: 
> kappa:=sqrt(K-l)/den: 



> Bl 

> B2 

> B3 



=kappa*(All-z*Al2+y*A13) 
=kappa*(A21-z*A22+y*A23) 
=kappa*(A31-z*A32+y*A33) 



In Maple, let us denote t 
norm 



as tell, tel2, tel3, and the Riemannian 



K - 1 



K 



- asB. 



We have: 

> tell:=u/alpha: 

> tel2:=v/alpha: 

> tel3:=w/alpha: 

> B:=sqrt(l-(1/K)): 

To reduce clutter, let us also introduce two abbreviations: 

> rho:=alpha/F: 

> phi:=(beta+alpha*B~2)/F: 

Now we are ready to give the formula for the inverse of the fundamental 
tensor. Denote that inverse, in Maple, as Gij. Then: 
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> Gll:=rho*All+rho"2*phi*tell*tell-rho'2*(tcll*Bl+tell*Bl): 

> G12:=rho*A12+rho~2*phi*tell*tel2"rho'2 s,! (tell s,! B2+tel2*Bl): 

> G13:=rho*A13+rho'2*phi*tell*tel3-rho^2*(tell*B3+tel3*Bl): 

> G21:=G12: 

> G22:^rho*A22+rho~2*phi*tel2*tel2-rho'2*(tel2*B2+tel2*B2): 

> G23:=rho*A23+rho~2*phi*tel2*tel3-rho'2*(tel2*B3+tel3*B2): 

> G31:=G13: 

> G32:=G23: 

> G33:=rho*A33+rho~2*phi*tel3*tel3-rho'2*(tel3*B3+tel3*B3): 



A. 4. The contravariant form of the geodesic spray coefficients. 

These are obtained by taking the covariant form Gi of the coefficients 
and raising the index with the inverse g lJ of the fundamental tensor. 

> Gl:=Gll*gl+G12*g2+G13*g3: 

> G2:=G21*gl+G22*g2+G23*g3: 

> G3:=G31*gl+G32*g2+G33*g3: 

A. 5. Getting set up for Berwald's formula. 

We assign names to the first and second order partial derivatives of G 1 . 
This will avoid having to re-compute them every time they are needed. Our 
Maple codes for Berwald's formula should run more efficiently as a result 
of this move. 

> Glx:=diff(Gl,x): 

> Glux:=diff(Glu,x): 

> Glvx:=diff(Glv,x): 

> Glwx:=diff(Glw,x): 

Duplicate the above with x replaced by y, z, u, v, w. 

Then duplicate all with Gl replaced successively by G2 and G3. 

A. 6. The spray curvature a la Berwald's formula. 

We now get Maple to calculate the nine components K l k of the spray 
curvature using Berwald's formula 

K \ = Z{G l ) x k - y J (G t ) x iy k - {G l ) y j{G : ')yk + 2G J (G l ) yJy k. 

Even though the covariant form Kit of the spray curvature is symmetric 
in its two indices, the same cannot usually be said of the type (J) form 
K l k :— g^Kjk- So, even at a purely numerical level, K\ ^ K l k in general, 
unless of course one is using a g-orthonormal frame. 

Let us use Kay[i, k] as the Maple names for those nine spray curvatures. 
The reason for not using K[i, k] is because of the next command. 

> Kay:=array(1..3,1..3): 
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Had we used K[i,k], the above would read "K:=array(1..3,1..3):" ; which 
would wreak havoc with our codes since the Maple variable K already 
stands for something else (namely the constant positive flag curvature). 

> Kay[l,l]:=2*Glx 

-Glu*Glu-Glv*G2u-Glw*G3u 

-u*Glux-v*Gluy-w*Gluz 

+2*Gl*Gluu+2*G2*Gluv+2*G3*Gluw: 

> Kay[l,2]:=2*Gly 

-Glu*Glv-Glv*G2v-Glw*G3v 

-u*Glvx-v*Glvy~w*Glvz 

+2*Gl*Glvu+2*G2*Glvv+2*G3*Glvw: 

> Kay[l,3]:=2*Glz 

Glu*Glw-Glv*G2w-Glw*G3w 
-u*Glwx-v*Glwy-w*Glwz 
+2*Gl*Glwu+2*G2*Glwv+2*G3*Glww: 

> Kay[2,l]:=2*G2x 

-G2u*Glu-G2v*G2u-G2w*G3u 

-u*G2ux-v*G2uy-w*G2uz 

+2*Gl*G2uu+2*G2*G2uv+2*G3*G2uw: 

> Kay[2,2]:=2*G2y 

-G2u*Glv-G2v*G2v-G2w*G3v 

-u*G2vx-v*G2vy-w*G2vz 

+2*Gl*G2vu+2*G2*G2vv+2*G3*G2vw: 

> Kay[2,3]:=2*G2z 

-G2u*Glw-G2v*G2w-G2w*G3w 

-u*G2wx-v*G2wy-w*G2wz 

+2*Gl*G2wu+2*G2*G2wv+2*G3*G2ww: 

> Kay[3,l]:=2*G3x 

-G3u*Glu-G3v*G2u-G3w*G3u 

-u*G3ux-v*G3uy-w*G3uz 

+2*Gl*G3uu+2*G2*G3uv+2*G3*G3uw: 

> Kay[3,2]:=2*G3y 

-G3u*Glv-G3v*G2v-G3w*G3v 

-u*G3vx-v*G3vy-w*G3vz 

+2*Gl*G3vu+2*G2*G3vv+2*G3*G3vw: 

> Kay[3,3]:=2*G3z 

-G3u*Glw-G3v*G2w-G3w*G3w 

-u*G3wx-v*G3wy-w*G3wz 

+2*Gl*G3wu+2*G2*G3wv+2*G3*G3ww: 



A. 7. The criterion for having constant flag curvature. 
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Finally, we ask Maple to check whether our Randers metric F has con- 
stant flag curvature K. The criterion we shall use has been derived in §6. 
It reads 



where r l k equals F 2 times the terms inside the parentheses. The pertinent 
Maple codes are: 

> Fu:=difT(F,u): 

> Fv:=diff(F,v): 

> Fw:=diff(F,w): 

> tau:=array(1..3,1..3): 

> tau[l,l]:=F~2-u*F*Fu: 

> tau[l,2]:= -u*F*Fv: 

> tau[l,3]:= -u*F*Fw: 

> tau[2,l]:= -v*F*Fu: 

> tau[2,2]:=F~2-v*F*Fv: 

> tau[2,3]:= -v*F*Fw: 

> tau[3,l]:= -w*F*Fu: 

> tau[3,2]:= -w*F*Fv: 

> tau[3,3]:=F~2-w*F*Fw: 

The condition we are striving for is K l k = Kr l k . To see if that's the case, we 
form the difference of the two sides, and also their ratio for good measure. 

> with(linalg): 

> f:= (i,k) -> Kay[i,k]-K*tau[i,k]: 

> dif:=matrix(3,3,f): 
Hopefully the answers are zero. 

> h:= (i,k) -> Kay[i,k]/(K*tau[i,k]): 

> quot:=matrix(3,3,h): 
These quotients should all be 1. 

A.8. The verdict. 

It was no problem for Maple to calculate symbolically those differences 
and quotients defined in §A.7. However, attempts to get Maple to symbol- 
ically simplify the answers to a or a 1 consistently crashed because the 
resulting expressions were too large. 

Given that, we did the next best thing. We randomly selected numerical 
values for the position coordinates x, y, z, the velocity variables u, v, w, 
and the positive flag curvature K. Then we asked Maple to evaluate those 
differences and quotients at the stipulated x, y, z, u, v, w, and K. In 
retrospect, it now appears that for the purpose of comparing K\ to Kr l k , 
forming the ratio of the two terms works better than taking their difference. 
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We noticed that the amount of RAM used in our computations is rou- 
tinely in excess of 1 Gigabyte. It is not clear whether this is attributable 
to any inefficiency in our codes. A sampling of our numerical results are as 
follows: 

> simplify(eval(dif[l,l],[x=1.0,y=2.0,z=-3.0, 

u=3.1416,v=2.78,w=137.0,K=29.0])); 

.009000000000 

> factor (simplify (eval(quot[l,l],[x=1.0,y=2.0,z=-3.0, 

u=3.1416,v=2.78,w=137.0,K=29.0]))); 

1.000000007 

> simplify(eval(dif[l,2],[x=9.0,y=7.0,z=-5.0, 

u=3.1416,v=2.78,w=137.0,K=31.0])); 

.008218000000 

> factor (simplify (eval(quot[l, 2], [x=9.0,y=7.0,z=-5.0, 

u=3.1416,v=2.78,w=137.0,K=31.0]))); 

.9999957122 

> simplify (eval(dif [1 ,3] , [x=9,y=7,z=-5, 

u=31416,v=278,w=137,K=31])); 



> factor (simplify (eval(quot[l, 3], [x=9,y=7,z=-5, 

u=31416,v=278,w=137,K=31]))); 

1 

> simplify(eval(dif[2,l],[x=131,y=17,z=-59, 

u=61413,v=872,w=l/137,K=2])); 



> factor (simplify (cval(quot [2, l],[x= 131, y= 17, z=-59, 

u=61413,v=872,w=l/137,K=2]))); 

1 

> simplify(eval(dif[2,2],[x=l,y=2,z=-3, 

u=71,v=5,w=l/137,K=29])); 



> factor (simplify (eval (quot [2 , 2] , [x= 1 ,y=2 ,z=-3 , 

u=71,v=5,w=l/137,K=29]))); 

1 

> simplify(eval(dif[2,3],[x=1.0,y=2.0,z=-3.0, 

u=3.1416,v=2.78,w=137.0,K=29.0])); 
-.005050000000 

> factor (simplify (eval(quot [2, 3], [x=1.0,y=2.0,z=-3.0, 

u=3.1416,v=2.78,w=137.0,K=29.0]))); 

1.000000215 
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> simplify (eval(dif [3 , 1] , [x= 1 ,y=2,z=3, 

u=5,v=7,w=ll,K=13])); 



> factor (simplify (eval(quot [3, l],[x=l,y=2,z=3, 

u=5,v=7,w=ll,K=13]))); 

1 

> simplify(eval(dif[3,2],[x=199.7,y=-2.4168,z=3.5, 

u=59,v=79,w=119,K=357])); 
.00003253970835 

> factor(simplify(eval(quot[3,2],[x=199.7,y=-2.4168,z=3.5, 

u=59,v=79,w=119,K=357]))); 

.9999963000 

> simplify (eval(dif [3, 3], [x=l,y=2,z=-3, 

u=71,v=5,w=l/137,K=29])); 



> factor (simplify (eval (quot [3 , 3] , [x= 1 ,y=2 ,z=-3 , 

u=71,v=5,w=l/137,K=29]))); 

1 

As one can see, the numerical evidence is overwhelmingly in favor of our 
Randers metric having constant positive flag curvature K. 



References 

[AIM] P. L. Antonclli, R. S. Ingarden, and M. Matsumoto (1993) The Theory of Sprays 
and Finsler Spaces with Applications in Physics and Biology, FTPH 58, Kluwcr 
Academic Publishers. 

[AZ] H. Akbar-Zadch (1988) Sur les espaces de Finsler a courbures sectionnelles con- 

stantcs, Acad. Roy. Belg. Bull. CI. Sci. (5), 74, 281-322. 
[Bl] L. Berwald (1947) Projektivkriimmung allgemeiner affiner Raume und Finslcr- 

schc Raumc skalarcr Kriimmung, Ann. Math. (2), 48, 755-781. 
[B2] L. Berwald (1941) Two-dimensional Finsler spaces with rectilinear extremals, 

Ann. Math. (2), 42, 84-112. 
[B3] L. Berwald (1926) Untersuchung dcr Kriimmung allgemeiner mctrischcr Raumc 

auf Grund des in ihncn hcrrschcndcn Parallclismus, Math. Z., 25, 40—73. 
[BCS] D. Bao, S.S. Chern, and Shcn, Z. (2000) An Introduction to Riemann-Finsler 

Geometry, Graduate Texts in Mathematics, vol. 200, Springer. 
[BM] S. Bacso and M. Matsumoto (1997) On Finsler spaces of Douglas type, a gener- 
alization of the notion of Berwald space, Publ. Math. Debrecen, 51, 385-406. 
[Brl] R. Bryant (1996) Finsler structures on the 2-sphere satisfying K = 1, Cont. 

Math., 196, 27-42. 

[Br2] R. Bryant (1997) Projectively flat Finsler 2-spheres of constant curvature, Se- 

lecta Mathematica, New Series, 3, 161—203. 
[Br3] R. Bryant, Finsler surfaces with prescribed curvature conditions, Aisenstadt 

Lectures, in preparation. 
[D] J. Douglas (1928) The general geometry of paths, Ann. Math. (2), 29, 143-168. 



FINSLER METRICS OF CONSTANT POSITIVE CURVATURE 



35 



[F] E.J. Flaherty (1976) Hermitian and Kahlerian Geometry in Relativity, Lecture 

Notes in Physics, vol. 46, Springer. 
[GLP] P. Gilkey, J. Leahy, and Jeong-Hyeong Park (1996) The spectral geometry of 

the Hopf fibration, J. Phys. A: Math. Gen., 29, 5645-5656. 
[M] M. Matsumoto (1989) Randers spaces of constant curvature, Rep. on Math. 

Phys., 28, 249-261. 

[R] H. Rund (1959) The Differential Geometry of Finsler Spaces, Springer- Vcrlag. 

[Ra] G. Randers (1941) On an asymmetric metric in the four-space of general rela- 
tivity, Phys. Rev., 59, 195-199. 

[Ro] C. Robles, private communication and doctoral dissertation in progress. 

[SI] Z. Shen (2000) Differential Geometry of Sprays and Finsler Spaces, IUPUI 
Monograph. 

[S2] Z. Shen (1996) Finsler manifolds of constant positive curvature, Cont. Math., 
196, 83-93. 

[YS] H. Yasuda and H. Shimada (1977) On Randers spaces of scalar curvature, Rep. 
on Math. Phys., 11, 347-360. 

Department of Mathematics, University of Houston, Houston, TX 77204- 
3476, USA 

Department of Mathematical Sciences, IUPUI, Indianapolis, IN 46202-3216, 
USA 

E-mail address: bao@math.uh.edu 
E-mail address: zshen<Smath.lupul.edu 



